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Abstract 

We propose a general formalism to study the static properties of a system com- 
posed of particles with nearest neighbor interactions that are located on the sites of 
a one-dimensional lattice confined by walls ( "confined Takahashi lattice gas" ) . Lin- 
ear recursion relations for generalized partition functions are derived, from which 
thermodynamic quantities, as well as density distributions and correlation functions 
of arbitrary order can be determined in the presence of an external potential. Ex- 
plicit results for density profiles and pair correlations near a wall are presented for 
various situations. As a special case of the Takahashi model we consider in par- 
ticular the hard rod lattice gas, for which a system of nonlinear coupled difference 
equations for the occupation probabilities has been presented previously by Rob- 
ledo and Varea. A solution of these equations is given in terms of the solution of 
a system of independent linear equations. Moreover, for zero external potential in 
the hard rod system we specify various central regions between the confining walls, 
where the occupation probabilities are constant and the correlation functions are 
translationally invariant in the canonical ensemble. In the grand canonical ensemble 
such regions do not exist. 

1 Introduction 

The understanding of the static and dynamic behavior of fluids in confined geometries 
is a problem of active current research [|l], ^. This research is largely motivated by 
technological applications where one wants to create small surface structures with suitable 
physical and chemical properties [0 . The question, how the formation of such structures 
is influenced by confining walls, has raised the interest in many basic phenomena, such 
as the development of density profiles, pair correlations and ordering effects at surfaces 
@], wetting transitions P], or a variety of surface induced kinetic processes 0. In 
general, an exact analytical treatment of the various effects is not possible and one has 
to rely on approximation schemes. Well established techniques for this purpose are the 
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density functional theory with its many variants (see ^ and references therein), the cluster 
variation |^ and path probability method |10|, as well as the classical thermodynamic 
perturbation theories (see e.g. ||TT[|). 

In one-dimensional fluid systems, however, exact results can be obtained for special 
models. Such results have their merits in establishing "boundary conditions" for the 
development of approximate theories in higher dimensions c? > 1 in the sense that these 
should become exact for = 1. Physically, this issue in particular pertains to theories for 
confined systems with tunable confinement, allowing a "dimensional crossover" form the 
case d > 1 to d = 1 [0. Exact findings also provide the possibility to systematically test 
the quality of approximations. Moreover, for many phenomena, such as e.g. the emergence 
of the well-known density oscillations of fluids near hard walls, the dimensionality does 
not seem to play a crucial role, and valuable insight into the origin of these phenomena 
may be gained by investigating appropriate one- dimensional reference systems. 

The first exact density functional in d = 1 was set up by Percus [jl3| for a fluid 
of hard rods. This functional yields an integral equation for the density profile in an 
arbitrary external potential. Later Percus showed that a further exact density functional 
can be written down for the special "sticky core" model iT^ that, in addition to the 



hard-rod repulsion, includes a "zero-range" attractive force between nearest neighbor 
rods. This was subsequently generalized to finite-range forces between neighboring rods 
| T5| . A generalized discrete version of the continuum hard- rod fluid on a linear chain 



was studied by Robledo and Varea ||I6|. They derived an exact functional for the mean 
occupation numbers of the rod centers on the chain, which, by taking the continuum 
limit, allowed them to recover the continuum density functional of Percus (for a review 
on classical density functionals, see Within this theory, the discrete hard- rod model 

leads to a rather complicated system of nonlinear finite difference equations for the mean 
occupation numbers in an arbitrary external potential, whose numerical solution requires 
a considerable computational effort. 

In this article we will show that a more general discrete one-dimensional system can 
be considered, which allows one to calculate bulk and surface thermodynamical proper- 
ties as well as equilibrium density profiles and density correlations in arbitrary external 
potentials. The system is an extension of the continuum Takahashi model |T^ to a lattice 
gas model, in which only neighboring particles interact with each other. We first show 
that the canonical and grand canonical partition functions of this "Takahashi lattice gas" 
(TLG) obey simple recursion relations and that density profiles and density correlations 
can be conveniently calculated from the partition functions due to the one-dimensional 
nature of the model. Since the TLG contains as a special case the hard-rod model stud- 
ied by Robledo and Varea |jl6| (for which the interaction potential would be infinite for 
interparticle distances smaller than the rod length and zero else), we can give a simple 
solution of the nonlinear difference equations derived in fl^. In this context we found it 



worthwhile to rederive the central formulae given in |jl6[. We will apply our formalism 
to a system of hard-rods confined both by hard and soft walls, and also to a system of 
particles, which in addition to the athermal hard-rod repulsion experience a finite inter- 
action potential over a limited range. For these different cases density profiles and density 
correlations near the confining walls will be discussed in detail and compared with each 
other. 
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2 Takahashi Lattice Gas 



In the TLG we consider N particles at positions i^, k = 1, . . . , N, on a linear chain with 
M sites i = 1, . . . , M. No more than one particle is allowed to occupy a given lattice site. 
Two neighboring particles separated by n — 1 vacant lattice sites interact via a potential 
v{n). There is no interaction between particles that are not nearest neighbors, that means 
between particles that have at least one other particle in between them. In addition the 
particles experience an external potential u{i) and it is assumed that two confining walls 
are present at the boundary sites i = and i = M + 1. These are modeled by two 
additional particles that are held fixed at the boundary sites. The energy of a particle 
configuration 1 < ii < . . . < ij^ < M is then given by 

N N 

PH = J2 u{ik) + v{ii) + J2 ^i^k - ik-i) + v{M + l-iN), (2.1) 

k=l k=2 

where (3 = l/ksT and we have assumed for simplicity that the particles at the boundary 
sites are the same as those on the chain. It will become clear in the following that one 
could also consider a modified interaction with the walls. 



2.1 General case of arbitrary external potential 

In the presence of a (non-constant) external potential u{i), it is convenient to define a 
generalized canonical partition function by 



Z{N, M', a)= exp - 

l<ii<...i]V<Af' 



N 



N 



^ u{ik + a)+ v{ii) + Yv{ik - ik-i) + v{M' + 1-%^) 

(2.2) 



.fc=i 



k=2 



for integers a > and M' + a < M. Equation (U) defines Z(iV, M', a) if < M', 
while for > M' we set Z(iV, M', a) = 0. Note that for a = and M' = M we recover 
the ordinary partition function Z{N,M) = Z{N,M,0). As far as only thermodynamic 
quantities shall be calculated, we could limit ourselves to the conventional form, but to 
evaluate density profiles and density correlations for the TLG, we need to consider the 
generalized functions (see below). Separating the summation over the positions I = of 
the rightmost particle, we can write 



M' 



Z{N,M',a) = ^exp [-f(M' + l -/)-«(/ + «)] X 

N-l 

exp - u{ik + a) + v{ii) + ^ v{ik - ik-i) + v{l - in-i) 

l<n<...<ijv_i<i-l k=l k=2 

M' 

^ exp [-f (M' + 1 - /) - m(/ + a)] Z{N - 1, / - 1, a) . (2.3) 



l=N 



N-l 



1=1 



In the last line we could start the summation from / = 1 because of our setting Z{N, M', a) i 
for A^ > M' . The recursion relation ( p.3|) is so far valid for A^ > 1. It becomes valid 
also for A^ = 1 if we set Z(0, M', a) = exp[-v{M' + 1)]. 
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From eq. ( |2.3| ) we readily derive a recursion relation for the analogous generalized 
grand partition function, 



Q{X,M',a) = ^Z(Ar,M',a)A^ (2.4) 

N=0 

M' 

= exp [-v{M' + 1)] + A ^ exp [-u{l + a) — v{M' + 1 — /)] Q{X, I — l,a) 



1=1 



As before, we define f2(A, M) = Q{X, M,0). This recursion relation (|2.5| ) is valid for 
M' > 1, but can be made valid also for M' = 1 if we set f2(A,0,Q;) = exp[— f(l)]. The 
fugacity A determines, for given M, the mean number A^(A, M) of particles in the system. 

Using the recursion relations ( |2.ij| , 23) one may calculate all thermodynamic properties 
for a given interaction v{n) and external potential u{i) by taking a = and M' = M. 
We now show how one can calculate also density profiles and correlations (of arbitrary 
order), once Z{N, M' ,a) or Q{X,M',a) has been calculated from eqs. ( p73| , p.5|) . To this 
end we first determine the probability w{l,r) to find the rth particle at position v = /, 



N 



N 



w{l, r) 



.<ir-l<l<ir + l< 



Z-\N, M) exp - v{ii) + v{ik - ik-i) +v{M + l-it,) + Y 



k=2 



k=l 



exp[-u{l)]Z-\N, M) 

r— 1 r— 1 

^(^i) + ^(^fc - ^fc-i) + ^^(/ - 1 + 1 - ir-l) + ^ u{ik) 



X ^ exp 

l<ii<...<V-i<i-l 



X Y exp 

l + l<ir + l<...<iN<M 



k=2 k=l 

N N 

v[ir+i-l)+ Y v{ik-ik-i) +v{M + l-iN) + Y "(^'^^ 

k=r+2 k=r+l 



exp[-u(/)] Z-\N, M) Z{r -!,/-!) 



N-r 



N-r 



X Y^ " 

l<jl<-<jN-r<M-l 



'(Jl) + Y ^(^''^ ~ •^''^-l) +HM-1 + 1- jN-r) + Y "(•^''^ + ^) 



k=2 



k=l 



exp[— u( 



Z{r -1,1- 1)Z{N -r,M -1,1) 
Z{N, M) 



(2.5) 



In the second step we have introduced the shifted particle positions js = ir+s — I for 
s = 1, . . . ,N — r. By doing this, ^(^a:) transforms to u{jk + ^^^1 it becomes clear now 
why we had to introduce the generalized partition function Z{N, M', a) with a ^ 0. In 
the canonical ensemble the probability p{l; N, M) for the site / to be occupied is then 
given by 

p{l; N,M) = f2 ^(^, r) = Z(r - 1, / - 1) Z{N -r,M- I, I) . (2.6) 



r=l 



r=l 
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In the grand canonical ensemble the convolution in eq. ( p^.6D factorizes and we obtain the 
corresponding occupation probability 

A, M) = A exp[-«(/)] ^ ^i^^) ' (2.7) 

We like to note that for a symmetric external potential, u[%) = u{M + 1 — i), it follows 
that n{X, M-l,l) = n{X, M-l). Hence it suffices to calculate ^(A, M', a) from eq. ( p3| ) 
for a = to obtain the density profile in this symmetric case. 

By an analogous decomposition of the partition function into products of generalized 
partition functions (corresponding to various system sizes) one can derive the joint prob- 
abihties Ps(Zi, . . . , iV, M) to find the sites li < . . . < Ig being occupied in the canonical 
ensemble, 

p.(/i,...,/.;iV,M) = ^^^^'^r^^'^^ 5^ {z(ri-l,/i-l)Z(iV-r„M-/„y 

^ ' ' l<ri<...<rs<N 



s-l 

X 

fc=l 

From this we obtain the corresponding joint probabilities in the grand canonical ensemble. 



s—^ 

JJZ(rfc+i - rfc - 1,/fc+i - Ik - 1,4)} • (2.8) 



p,(/i,...,/,;A,M) = ''''\fh) ""^^'^^ ~ ~ 



J]fi(A,4+i-4-l,4). (2.9) 



-1 

X 
A:=l 

Note that for s = 1 eqs. (|]8|, U) reduce to eqs. From eqs. (U, ^[SD or eqs. (U, 

|2.8| ) one can readily calculate density profiles and density correlations of arbitrary order 
in the canonical or grand canonical ensemble for arbitrary interaction v{n) and external 
potential u{i). 

2.2 Special case of vanishing external potential 

In case of a vanishing (or constant) external potential it is not needed to introduce the 
generalized partition functions in eqs. (p.2| , |2.5| ) and accordingly one can set the third 
argument a in f2(A, M, a) equal to zero in all formulae in Sec. |2.1| . The occupation 
probabilities p(/; A^, M) and p{l] A, M) in eqs. ( |2.6| , p.7|) can be written as 

, A , , Z(N-r,M-l) , n(X,l-l)n(X,M -I) 
p(l; iV. M) = g Z(r - M - 1) ' . Pd: A. M) = ^^^-J^^ . 

(2.10) 

and analogous simplifications are obtained for the joint probabilities of higher order in 
eq. (^). 

Moreover, we can solve the recursion relations ( p.3| , |2.5|) explicitly in terms of the 
generating functions H{N,s) = ^f^^^ Z{N, M)s^ and G{X,s) = Em=o ^(^' ^)^*^' 
which are explicitly given by 

H(N,s) = ^^. G(A..) = -4' (2.11) 

s s[l — X(p[s)\ 
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with 



ip{s) = J2^xp[-v{l)]s'. (2.12) 



1=1 



If v{l) has a finite range, that means v{l) = for I > Iq, we obtain from (|2.11|) ^(A, s 



P(A, s)/Q(A, s), where P(A, s) and Q{X,s) are polynomials in s of degree Iq — 1 and Zq; 
respectively. According to a theorem for rational generating functions pO[, f2(A, M) then 



has the form fi(A, M) = ^J^g ^(A, M)sj{X)~^ , where Sj(A), j = 0, . . . , A; are the distinct 
zeros of Q{X,s) with multiplicities dj, and Cj(A,M) are polynomials in M of degree less 
than dj. The moduli of the zeros Sj are considered to be ordered, |so| < < ... < \sk\- 
As shown in Appendix A, sq is real with < sq < 1, cio = 1, and \sj\ > sq for 
j = 1, . . . ,k. Hence we can write i7(A, M) = Sq^'^[cq + X]j=i Cj(^)(so/sj)*^] and obtain 

fi(A,M) ~co(A)so(A)-*^ (2.13) 

in the thermodynamic limit M —>■ oo. The one-to-one correspondence between the fugac- 
ity A and the number density p = N{X, M)/M in this limit follows from the relations (see 
Corollary A.l. in Appendix A) 

p = ^{so)/[so^\so)], X = l/^{so). (2.14) 

Using the asymptotic limit for f2(A, M) we obtain from eq. ( p.lO| ) 

p^{l; A) = lim p{l; A, M) = A^](A, / - 1)4 . (2.15) 

M^oo 



In fact, as shown in Appendix A, eqs. ( ^.14| , |2.15| ) hold true even for a more general 



interaction potential v{l), which for I larger than some is bounded and for / oo 
approaches zero. The analogous occupation probability Poo(0 canonical ensemble 

is the same as Pooi^', A), if for given p the corresponding unique fugacity A is used (see 
eq. ( |2.14|) and Appendix A). Moreover, the joint probabilities Ps{h, • • • A, M) in the 



grand-canonical ensemble (and the corresponding ps{li,... ,ls',N,M) in the canonical 
ensemble) factorize in terms of Poo{l', A) in the thermodynamic limit, that means 

s 

Ps,oo{h, ■■■ ,h;X) = lim ps{h, . . . ,ls;X,M) = Poo{h\ A) TTj9oo(^fe - h-i, A) . (2.16) 

k=2 

For the special case of a finite range interaction potential considered above (i.e. v{l) = 
for / > Zq) there exists a constant C > such that \paoil', X) — p\ < Cl'^e~''^^, where 
< Zo — 2 is an integer, and ^ = — 1/ ln(r) with r = ma.xi< j<k{so/\sj\} < 1. 

3 Hard-Rod Lattice Gas Revisited 

A particularly simple situation occurs, when the interaction potential in eq. ( p.l| ) is given 

by 

oo , < n < 2m 
v{n) = fHR(n) = { (3.1) 
, n > 2m 
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This potential can be viewed as describing a system composed of hard-rods with lengths 
2m (with hard walls at positions m and M + 1 — m due to the fixed hard- rods at positions 
and M+ 1). 

Setting the mass density of the rods equal to unity, we can express the local mass 
density Pmass(^; M) along the one-dimensional chain by the occupation probabilities 
N, M) (that refer to the rod centers) according to 

^ m—l 

p^Uh M) = ^W- M) + p{l + m; iV, M)] + ^ p(/ - j- N, M) . (3.2) 

j=-(m-l) 

This formula holds true in the canonical as well as in the grand canonical ensemble (if 
Pmass(^; M) is replaced by Pmass(^; A, M) and p{l] N, M) by p{l] A, M)). Note that the 
total number N of rods must be smaller than M/2m. 



3.1 Explicit results for homogeneous systems 

Density profiles and density correlations can be calculated explicitly in the absence of an 
external potential by using the general method developed in Sec. |2]^. From eq. (p.ll|) we 
find ^{s) = s2"/(l - s) and G(A, s) = s^'^'^^/il - s- As^'"), and therefore 

n(A. D^ti'- <^"' - ^ A" . ZiN, ,)J'- (^'" -;)(^ + . (3.3) 

n=0 V / \ / 

The occupation probabilities p{l] N, M) and p{l; A, M) in the canonical and grand-canonical 
ensemble then follow by inserting these expressions into eqs. (|2.6| , p.lOD , and the correla- 



tions analogously. One can show |jT8| that p{l; N, M) and p{l] A, M) become maximal at 
the points / = 2m and / = M + 1 — 2m closest to the walls. The reason for this is that 
by fixing the position of a rod next to a wall the number of possible configurations (and 
hence the entropy) for the remaining (A^ — 1) rods will be largest. 

Extending this line of thinking one would guess that the most likely configuration near 
a wall is that where the rods are at positions / = 2m, 4m, 6m, .... One then should expect 
oscillations in the occupation probabilities to emerge with a period of typical size 2m. In 
fact, in the thermodynamic limit M — > oo one finds A = 1/'^{sq) = (1 — so)sq^™ with 
So = (1 — 2mp)/[l — (2m — l)p\ (see eq.( p.l4|) ). Moreover, as shown in Appendix B, the 



zeros Sj = \sj\ exp{i6j) (see Sec. p.2| ) are all different, and eq. ( |2.15| ) becomes 

2m— 1 X \ 

i=i ^i^^i^ 



2m-l y X I 

Poo(/,A) =P+T c,(A) ^ e-^^'. (3.4) 



The dj are in the open interval < 6j < 2tt, that means the profile Poo(^ A) is a super- 
position of simple oscillating and exponentially decaying functions. When considering a 
system of finite length, the effects induced by the second wall at position M + 1 — m 
on the profile near the first wall at position m are of order l/M in the grand-canonical 
ensemble and of order /^/M in the canonical ensemble. This is proven in Appendix B 
(more precisely, / must be of order o(M^/^) in the canonical and of order o(M) in the 
grand-canonical ensemble to obtain vanishing contributions in the thermodynamic limit.) 
Hence, the finite size corrections to eq. (p.4|) become small for large M. 
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More surprising, if the number density is smaller than half of that for the closed 
packed configuration, i.e. p < l/4m, one can show (see Appendix B) that in the canonical 
ensemble the p{l; N, M) are constant (/-independent) inside the central region IZi = {h E 
N I < k < M + 1 - /(I)} with /(I) = (2m - l)N + 1. At the outer boundary 
points = - 1 and = M + 2 - p{P; N, M) is different from the constant 
value inside T^i, that means TZi is "maximal" in the sense that there exits no other 
region of constant occupation probability enclosing parts of TZi. It is interesting to note 
that |p(/(^) ± 1;A^,M) - p{l^^^ ; N , M)\ = l/Z{N,M), that means the logarithm of the 
jump in the occupation probability at the boundaries of TZi provides the free energy oc 
log2'(A^, M). Furthermore, the joint probabilities Ps{h, ■ ■ ■ N, M) are translationally 
invariant inside ("maximal") regions TZg = {{h, . . . ,ls) E W) \ < /i;2m < h — 
h-i for = 2, . . . , s; < M + 1 - with = (2m - 1)(A^ + 1 - s) + 1,0 that 
means there exists a function /(xi, . . . , Xg^i, N, M) such that for all (/i, . . . , /g) G IZg, 
Ps{h, ■ ■ ■ N, M) = f{l2 — li, . . . ,ls — h-i] N, M). Corresponding regions have been 
found in the continuum version of the hard-rod lattice gas, the so-called hard-core fluid 



model |2T|. Moreover, if (Ij+i-lj) > [(2m- 1)(A^-M -s) -t- 1] for all j = 1, . . . then 
Ps{h, ■ ■ ■ Js', N, M) is constant for (/i, . . . Js) £ T^s (this was shown to hold true in the 
continuum model for the pair distribution functions only [Q). As is shown in Appendix 
B also, the situation is quite different in the grand-canonical ensemble. Here, there exist 
no regions of constant occupation probabilities and translational invariance of the joint 
probabilities (except for trivial cases^). 



3.2 Free Energy Functional 



An alternative way to treat the hard-rod lattice gas has been followed in [16]. In this 
approach, which relates to density functional theory of classical fluids, one considers 
the grand-canonical ensemble and defines the occupation numbers Xi, where Xj = 1 if 
site i is occupied by a rod center, and Xj = else.0 Note that these random vari- 
ables are not independent: Since the rods have size 2m we have to require Xj = 
for \j — i\ < 2m if Xj = 1. We define Cm as the set of all allowed configurations 
{xj}. The idea then is to calculate explicitly the probability xi^iy • • • ; ^m) of an allowed 
configuration (xi, . . . ,xm) ^ Cm by regarding the occupation probabilities pi = (xj) 
as fixed (([...]) = ,XM)eCAf [• • • • • • '^m))- For given pi it turns out that 

logx(a;i, . . . ,xm) depends linearly on the occupation numbers (which is a fortunate fea- 
ture of the hard rod system, see below). By equating xi^i^ • • • ? ^m) with the Boltzmann 
formula for all (xi, . . . , xm) £ Cm we have 



M 



logfi(A,M) = logx(a;i, . . . ,xm) + ^Xs[u{s) - fx] , (3.5) 



^By saying that the Ps{h, ■ ■ ■ ,ls;N,M) are not translationally invariant if /i < . . . < ^ 7?.s we 
exclude the trivial case, where h and Is must be occupied by the first and last rod center, respectively 
(i.e. for li < 4m and Is > M + 1 — 4m). (In this case, our system can be considered as being composed 
of — 2 rods on a chain of length {Is — h) and the ps{li, . . . ,ls',N, M) are translationally invariant then 
on trivial reasons.) 

^For the following it is convenient to make the transformation I ^ I — (2m — 1) of the site positions 
and to change the system size according to M ^ M + 4m — 2. After these replacements the rods at the 
boundaries are at positions —2m + 1 and AI + 2m corresponding to hard walls at positions — m + 1 and 
M + m. Accordingly, the possible positions of the rod centers are 1, . . . , M. 
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where fi = log A is the chemical potentiaL Taking now the expectation value of ( ^.Sj ) with 
respect to the Xi, an exact density functional pj-'{pi, . . . ,Pm) = — logf2(A,M) of the pi 
is obtained. 

In order to determine • • • j^m) for (xi, . . . ,xm) G Cm we will make use of a 

Markov property. (Note that the Markov property is not valid in the canonical ensemble.) 
The constraint given by the finite rod lengths implies that the conditional probabilities 
Ws{xs\xs_i, . . . Xi) for the occupation number at site s to be Xg, given Xg^i, . . . Xi, are 
independent of Xi, . . . Xs_2m- In other words Ws{xs\xs-i, ■ ■ ■ Xi) fulfills the generalized 
Markov condition (which is rather obvious here but can be proven rigorously too [Q) 

Ws{Xs\Xs-l, . ..Xi) = Ws{Xs\Xs-l, . . . , Xs-2m+l) ■ (3.6) 

Due to this property we can express the joint probabilities xi^i^ ■ ■ ■ ^^m) as 

X{Xl, . . . ,Xm) = Wi{Xi)w2{x2\Xi) . . .W2ni{x2m\x2m-1, ■ ■ ■ ,Xi) ... 

XWs{Xs\Xs-l, ■ . . ,Xs-2m+l)---'WM{xM\xM-l,--- ,XM-2m+l ) (3.7) 

It is convenient to formally extend the system to integers i < and to set Xj = for all 
—2m + 2 < i < 0, such that we can write eq. ( |3.7|) in the compact form 

M 

XiXi, ... , Xm) = Y\_ Ws{Xs\Xs-l, . . . , Xs-2m+l) ■ (3.8) 
s=l 

For calculating Ws{xs\xs-i, . . . , Xs-2m+i) we have to deal with two cases only: (i) One 
of the given random variables Xs-2m+i, ■ ■ ■ , Xg-i is equal to one and the rest of them equal 
to zero, and (ii) all Xs-2m+i, ■ ■ ■ , Xg^i are zero. In all other cases there would be at least 
two of the Xs-2m+i, ■ ■ ■ , Xs-i equal to one, but this is not allowed, because it would imply 
that rods overlap. For the same reason we must have = in situation (i), that means 
we obtain 

Ws{Xs = 0\Xs-l = 0, . . . ,Xi = 1, . . . , Xs-2m+l = 0) = 1 . (3.9) 

The situation (ii) is more complicated. By definition we can write 

wsixs\0, . . . , 0) = — -— , (3.10) 

l^s-l,s-2m+l\Xs-l — U, . . . ,Xs-2m+l — UJ 

where ni^ki^i, . . . ,Xk) is the joint probability for the configuration {xi, . . . ,Xk} to occur. 
For k — I < 2m the normalization condition yields (again because of the non-overlapping 
condition) 

I 

1 = ^ Ki^k{xi, ... ,Xk) = K;,fe(0, ... , 0) + ^ /t/,fc(0, ... , Xj = 1, . . . , 0) . (3.11) 

{xj} j=k 

By definition we further have (for k — I < 2m) 

Pi = ^Xi Ki^k{xi, ... ,Xk) = fi:«,fc(0, . . . , Xj = 1, . . . , 0) , (3.12) 
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and hence it follows form eqs. ( p.ll| , p.l2D , 



Inserting the i^i^kixi, 



K/,fe(0, ... ,0) = 1 - . 

j=k 

,Xk) from eqs. ( |3.12| , |3.13| ) into eq. ( |3.1CI| ), we obtain 
,0) = 



(3.13) 



wJx.M, 



i-tm{s) 

i-C(^) ' 

Ps 



0: 
1, 



(3.14) 



where we have defined tmis) 



E 



2m- 1 ~ 

i=o 



Ps^j and C(s) 



E 



results ( |0| , p.l4| ) can be combined to express w(a;s|a;s_i. 

[1 tmi^i 



2m— 1 ~ i , , 

j=l Ps-i = ^m[S) - p. 



= tm{s) - Ps- The 
Xs-2m+i ) in the general form 



w{x g\x g—\^ ... , X 



s-2m+l 



(3.15) 



Note that the Xi appear linearly in the exponents of the transition matrix (p.l5 ), such 
that by inserting (|3.15|) into eq. ( p.8|) , and taking the logarithm, we find that logx is 
linear in the Xi. Using eq. (|3.5|) and averaging over Xi we finally obtain ||16 



M 



M 



(3T{pi, . . . ,pm) 



^ Ps [u{s) - ^i\ + ^Ps log p. + 

s=l s=l 
M M 

J](l-t„(s))log(l-t™(s))-5^ 



(3.16) 

:i-c(s))iog(i-c(s)) . 



s=l 



The functional ( p. 17 ) becomes minimal for the equilibrium density profile pi = p{l; X, M). 
The corresponding system of equations reads (/ = 1, . . . , M) 



^^(pi 
dpi 



,Pm) 



-l2 + u{l)+\0gpi+ J2 l0g(l-L(5))- 



s=l+l 



l+2m-l 

- E 

s=l 



log(l - t™(s)) = 
(3.17) 



It is clear that p{l; A, M) from eq. (|2.7|) (after making the transformations I I — {2m — 
1) and M — »• M + 2 — 4m) must solve ( p.l7| ). From a mathematical point of view 
this is an interesting example, where a system of coupled nonlinear difference equations 
(eq. ( p.l7| )) can be mapped by a nonlinear transformation (eq. ( p.7| )) onto a simple system 
of independent linear difference equations (eq. (|2.5| )). For the special case of a vanishing 
external potential even an explicit solution exists (see eqs. ( [^.10| , |3.3D). A direct proof 
that p{l\ A, M) from (|2.7|) indeed solves ( 3.17| ) is given in Appendix C. 

Next we rederive the exact free energy functional of Percus |jl3| by taking the proper 
continuum hmit of eq. ( p.l7|) . To do this we first have to note that eq. (|3.17| ) gets 
slightly modified, when it is viewed as resulting from a discretized form of an originally 
continuous system. This continuous system is defined by hard rods of length a with 
positions < yi < L, T/j+i — yi> a. In a discretization, we may subdivide the continuous 
system into M intervals X, {s = 1, . . . , M) of equal size Ay = L/M, and may set the 
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rod length 2m in the new discrete variables s equal to the integer part of aL/M. The 
occupation number Xg of the interval Xg is defined to be zero, if none of the yi G Tg and 
one else. Then the joint probability q{ii, . . . ,iN) to find rods at positions yi, . 



in the intervals Tj^ , . . . , Xj^ is, 



k = l -^^ik 

/ N ^ 

fi(A, L)-\Ayf exp - J2[u{ik) - /x] 



fc=i 



Since there is a one-to-one correspondence between the sets {ii, 
we immediately obtain 



Ay 
zat} and {xi 



M 



(3.18) 
, • • • ,xm} 

(3.19) 

Repeating the steps leading to eq. ( p.l7|) we get a modified functional J^{pi, . . . ,Pm), 
which is the same as given in ( p.l7|) plus the term [—{^^ps)\ogAy + o{Ay)/ Ay]. The 
ps are related to the occupation number density p{y) in the continuous system by Ps = 
jj^ dypiy) = [p{sL / M) Ay + o{Ay)] and by inserting this in the modified form of eq. (|3.17|) 
we obtain in the limit M — * oo {Ay —>■ 0) the Percus functional 



X{xi,... ,xm)=^{KL) ^(Ay)^-i"^exp -^[n(s) -/i], 



s=l 



1 + 



o(Ay) 
Ay 



dyp{y) My) - p + \ogp{y) - [1 + log(l - t{y))]} 



(3.20) 



where t{y) = f^^^dz p{z). To consider eq. ( |3.20|) as a mass density functional one should 

remember the relation pmass(l/) = t{y + o"/2) = jy^^j^ dzp{z) between the mass and the 
number density (that might be easily inverted by Laplace transformation). 



4 Density Profiles and Pair Correlations Near Walls 

In this Section we calculate density profiles and correlations for some cases to exemplify 
the formalism developed in the previous Sections ^ and |^. 

Figure 1 shows the occupation probability Pooi}) for a system of hard rods as a function 
of the distance / from a hard wall for (a) p = 0.1 and various (half) rod lengths m = 2, 3, 4, 
and (b) p = 0.02 and m = 14, 18, and 22 (poo(0 '^^s calculated from eqs. (p3|J2.15| )). As 



can be seen from the figure, Poo{l) exhibits oscillations with a period of order 2m, which 
become more pronounced with increasing m. For large /, Poo(0 approaches p. Note that 
for p = 0.1 the closed packed situation occurs already at m = 5 and the discreteness 
of the system is important (see Fig. la), while p = 0.02 (Fig. lb) corresponds to a 
continuum situation. The data in Fig. lb indicate that Poo{l) in the continuum limit (see 
Sect. ^.21 ) might have a discontinuity in the first derivative at the first minimum. Indeed 
this discontinuity occurs and its origin can be understood from the solution of the discrete 
system: From eq. (|2.15|) and the recursion relation (P3|) one derives 



Poo(0 = soPoo{l - 1) + (1 - sq)poo{1 - 2m) . (4.1) 
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Accordingly, when / < 4m, the second term in ( |4.1| ) is zero up to the first minimum in 
Poo(0 / = 4m — 1, and it first contributes when / = 4m. The additional contribution 
from the second term yields the discontinuity in the first derivative. 
The correlation function 

C{1) ^ P2,oo(2m, /) - poo(2m)poo(0 (4.2) 

between the first possible position 2m of a rod center and another rod center that is at 
distance / from the wall is shown in Fig. 2 for the same parameters as in Fig. 1. Similar 
as Poo{l), C{1) oscillates as a function of / with a period of order 2m; the strength of the 
oscillations increases with increasing m. For large Z, the absolute values of C{1) at its 
local maxima and minima decrease exponentially with /. 

Next we calculate density profiles and correlations for more general cases. To this end 
we consider (i) hard rods (f (/) = fHR(O) the presence of a "soft wall" with an attractive 
potential 

(5ua{l) = -5exp(-//20) , (4.3) 
and (ii) particles with a Lennard- Jones type Takahashi interaction of the form 

oo / < 2m 

pv^j{l) = { -4 2m < / < 3m (4.4) 
else 

in the presence of a hard wall (m(/) = 0). Figure 3 shows Poo(0 these two cases in 
comparison with the hard rod system for (a) m = 4 and (h) m = 18 (to calculate Poo(0 
for u{l) = uq{1) we have chosen a large system size M = 10^ and used eqs. ( p.5| , p.7|) . For 



both cases oscillations occur similar as in the hard rod system. In Fig. 3a the attractive 
wall potential causes the maxima and minima to become more pronounced than in the 
other cases, while in Fig. 3b only the occupation probability for the first rod next to the 
wall is strongly enhanced. Because the first rod center is strongly attracted by the wall, 
the position of the following minima and maxima of poo are shifted toward the wall. The 
weaker effects of the external potential in the continuum-like situation (Fig. 3b) are due 
to the fact that the first minimum of Poo(0 occurs at a position, where Uo{l) is already 
very small. For the Lennard- Jones type interaction we find the oscillations in Fig. 3a to 
be stronger than in the hard rod system, but the probability of the first rod to be right at 
the wall is reduced (for smaller m, however, p^{2m) can be larger than in the hard rod 
system). As can be seen from Figs. 4a,b, the changes of the correlation functions caused 
by the external potential uq{1) and by the interaction potential fLj(0 are fully analogous 
to the changes found for Poo(0 Figs. 3a,b. 



5 Summary 

As demonstrated in Sec. the recursion relations derived in Sec. ^ provide an efficient 
method to determine density distributions and correlations in the Takahashi lattice gas 
for arbitrary interactions and external potentials. We have proven in Appendix A that 
the wall-induced density oscillations decay exponentially into the bulk, if the interaction 
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potential has a finite range. If tlie potential has no finite range but still decays to zero 
when the inter-particle distance increases toward infinity, one can show only that the 
density for large distances from the wall converges to a constant bulk value. 

In the special case of the hard rod lattice gas, we have rederived the exact free energy 
functional of the occupation probabilities and the associated nonlinear system of coupled 
nonlinear difference equations. The general formalism derived in Sec. |^ allowed us to give 
a solution of these nonlinear difference equations in terms of the solution of a system of 
independent linear equations. Furthermore, various central regions between the confining 
walls have been specified in the canonical ensemble, where the occupation probabilities 
are constant, and where the correlations functions are translationally invariant or even 
constant too. In the grand canonical ensemble it was shown that such regions do not exist 
(except for trivial cases, see above). 

It is possible to extend the calculations for the hard rod system to more general 
situations, as, for example, to systems where certain groups of particles have differing rod 
lengths, or even to systems with randomly distributed rod lengths. 
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APPENDIX A 

For |s| < 1 and A G (0, oo) the generating functions of the canonical and grand- 
canonical partition functions are (compare with eq. ( |2.11|) ) 



(A.l) 

with 

oo 

^(s) = ^exp[-i;(/)]s', (A.2) 



1=1 



where v{l) has the following properties 



(i) max \v(l)\ < oo for some < oo , 

^ ^ h<l<oo 

(ii) lim v{l) = . 



(A.3) 



Lemma A.l. (See also |^.) For given A G (0, oo), there exists exactly one real 
positive solution so(A) G (0, 1) of 1 — Xf{s) = and this root is simple. If Sj G C is 
another root of 1 — Xip{s) = 0, then \sj\ > sq- 

Proof. The series Yl'ili I exp[— t>(/)]s'| converges for |s| < 1, and accordingly (p'{s) = 
^^j^ Z exp[— t>(/)]s'~^ > for positive real s. Since ^^(O) = and ^^(l) = oo, there 
exists exactly one sq G (0,1), for which (p{so) = 1/A (with multiplicity one). If sj 
were another root of 1 — Xip{s) = with \sj\ < sq, then |</3(sj)| < v'(lsjl) < v^(so) = 
1/A = (p{sj), which is impossible. If Sj = soexp(i^) with 9 G (0,27r) we have to require 
X^i^i 6xp[— f (/)]sq[1 — cos(/^^)] = and hence 6 = 27ck, A; G Z, in contradiction to the 
assumption 9 G (0,27r). 

Theorem A.l. Let So(A) G (0, 1) be the unique real positive solution of I — Xip{s) = 
0. Then 

fi(A,M) ~ co(A)so(A)-*^, for M-^oo. (A.4) 

where co(A) = ip{so)^ / slip'{so). 

Proof. Let ak = Q{X,k — 1)sq and fk = exp[— t'(A;)]sQ/ ^^^-^ exp[— f (/)]sq, for k = 
1,2,..., and set A = l/ip{so). Then, according to eq. ([^.5[), 

M-l M 

aM = V2(so)/m + ^ fu-iai = ^ au-ifi , (A. 5) 

1=1 1=1 

where Oq = <^{so). Obviously, fi>0 and Yl'^i fi ~ ^- employing the renewal theorem 



25| it follows 



«Af ~ ^<^° jr for M ^ OO . (A. 6) 

1^1=1 

Since YaLJIi = SoV5'(so)/</'(so), we obtain f](A, M) = aAf+iSg ^^^"^^^ ~ v5(so)^So ^'^/sg(/)'(so) 
Co(A)so for M ^ oo. 

Corollary A.l. 
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(i) poo{l; A) = \imM^ooP{l; A,M) = XQ{X,l - 

(ii) Ps,oo{h,--- ,^s;A) = limAf-^ooPsl^i, • • • Js;\M) = Poo{h] \) 111=2 Podh - h~u 

(Hi) p = (p{so)/ SQip'^so), where p = limM-^oo N{X, M)/M. {N{X,M) is the mean number 
of particles for given fugacity A and M.) 

Proof. According to eq. (CT) , p{l; A, M) = Afi(A, / - 1)^](A, M - /)/^](A, M) such 
that p{l; A, M) ~ Afi(A, / - l)co(A)so ^*'^~'Vco(A)so *^ = A^](A, / - 1)4 by theorem A.l. In 
particular, Poo(^;A) = \Q{X,l — 1)sq. Analogously, using eq. (|2.9|) , Ps,oo(^i7--- Js', ^) = 
A'*il(A, Zi - 1) Yll=2 ^(A, 4 - 4-1 - 1)4" • Since 4 = /i + Y.l=2i^k - h-i), we can write 

Ps,ooih,--- ,4; A) = Al](A,/i - 1)4 nL2[^^(^'4 - 4-1 - 1)4"''""'], which together 
with (i) gives (ii). By definition and theorem A.l, N{X,M) = XdlogQ{X, M)/dX ~ 
A91og[co(A)so(A)"^^]/9A, from which follows p = — Aso(A)/so(A). But from Av9(so(A)) = 1 
we immediately obtain — Aso(A)/so(A) = ip{so)/ soip'{sQ) and hence (iii). 

Theorem A. 2. Let Nm be any sequence with limM^oo Nm/M = p. Then 
(i) p^{l) = \imM^ooP{l;NM,M) =Poo(^; A) = Afi(A,/ - l)s|) 

(li) limM^ooPsih,- ■ ■ Js]Nm,M) = Ps,ooih, ■ ■ ■ ,4; A) = Poo(4; A) nL2Poo(4-4-i; A) , 

where A = l/ip{so) and sq is the unique positive solution of ip{so)/soip'{so) = p. 
Proof. In order to derive the asymptotic limit of the occupation probability 

p(l-,N.uM)^±Z(r-lJ-l) '^%'-'^'-'^ (A.7) 

in the thermodynamic limit, we use (see theorem 6.1 in 



ZiNM-l,M-l) = M?''\ u/2 [^ + Q(^'')] M^oo, (A.8) 

where = dl[\ogip{sMe'^)]u=o and sm is the unique non-negative solution of Nm/M = 
<f{s)/sip'{s). 

We further define sm as the unique non-negative root of {Nm—t) / (M—l) = ip{s) / s{p'{s) 
(for r,l given integers) and af^ = dl[\ogip{sMe^)]u=o- Then it is easy to show that there 
exist sequences Pm and 7m converging to finite values for M ^ oo with the property 



Sm — Sm 



1 + ^ + 0(M-2) 



~ 

<^M — 



l + ^ + 0(M-2) . (A.9) 



(For example, [3m = i} — r)f{sM)/sMf'{sM) with /(s) = ip{s)/sip'{s) has the desired 
properties.) Replacing Nm by Nm — r, M hj M — I, sm by sm, as well as ctm by ctm in 
eq. ([A.8|) , and using eq. (|A.9|) , we obtain after simple calculations 

Z{Nm -r -1,M -I -1) = Z{Nm -1,M ~l)s[MsMy"[l + 0{M-^)] for M ^ oo . 

(A.IO) 

Taking the limit M oo in eq. (|A.7|) we thus get (note that limM^oo sm = sq) 

I 

lim p{l; Nm,M) = y" Z{r -1,1- l)(^(so)""4 = Af^(A, / - 1)4 . (A.ll) 

r=l 
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The proposition (ii) follows by using the asymptotic form of Z{Nm —r,M — I) (eq. ( |A.10| ) 
in formula (|2.8|) ). 



APPENDIX B 

For |s| < 1 and A G (0, oo) let 

P(\ q\ °° „2m-l 

be the generating function from eq. ( [^.11| ) for the special case of the hard rod interaction 



potential defined in eq. (|3.1| ). 

Lemma B.l. The roots Sq, . . . , S2m-i of the polynomial Q{X, s) are all distinct. 

Proof. For one of the roots Si not to be simple, we must require that both Q{X, Si) = 
and idQ{X,s)/ds)s=s, = 0. But if idQ{X, s)/ds)s=s, = 0, we have s^^^-^ = -l/2mX and 
inserting this result into Q{X,Si) = we obtain Si = 2m/ {2m — 1) > 0, i.e. a positive 
real number. On the other hand, the only real solution of sf^~^ = —l/2mX is negative, 
which is a contradiction. Hence all roots must be simple. 



Let 

2m-l . V I [l/2m] . x s 



/I \r l—2mr 

(1 - So) So 



be the occupancy probability in the thermodynamic limit, where [x] denotes the integer 
part of X (see eq. (P^, and eqs. ( |2.14j , |2.15| , PT^D , and note that (p{so) = Sq™/(1 — sq)). 

Lemma B.2. For A^o = l/<^(so) = (l-So)s~^"' with sq = (1 -2mp)/[l - (2m- 
p = N/M {0 <p< l/2m), and P/M^ = o{M~^) 

(i) p{l- N, M) = p^{l; Aoo)[l + 0{iyM)] 

(ii) p,{h, ... ,h = l;N,M)= p^^^ih, ...,/, = /; Aoo)[l + 0(/VM)] 
Proof. According to eqs. (|2.(j|), 

p(/;iV.M)^|:z(.-M-l) ^'^^-;;;;-" . (B.2) 

With the definition 

where u = 1 — (2m — l)p, and x, y are integers, we can write (see the results for Z{N, M) 
in eq. 



^^^Z{N,M) = "''P ^^2"" - - 1) - I, -r) - f {-{2m - 1), 0)] . (B.4) 
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If x,y = 0(1) we obtain, by applying Stirling's formula, n! = (27rn)-'^/^(n/e)"exp[0(l/n 
f{x,y) = ^log( „ , ^ ) + M[u\ogu - {u -p)\og{u -p) -plogp] 



2 ^\2Tip{u-p)M 



+ xlogu - {x - y) \og{u - p) - y logp + O ( — ) (B.5) 



Note that the sum over r in eq. ( p.2|) runs at most up to the integer part of //2m, such that 
the arguments of the / functions appearing in eq. ( [B.4|) are all of order 0(/). Accordingly, 



Z{N -r,M -I) 
Z{N,M) 



;1_,„)-4^-2-)[i + o(/VM)], (B.6) 



from which we obtain (i) hj using eq. ( [B.2| ). 
Analogously, starting with eq. ( p.8|) . 



Psih,...,ls;N,M)= J2 ^(ri-l,/i-l)n^(r,+i-r,-l,/,+i-/,-l) ^^^J^'^ , 

l<ri<...<rs<N k=l ^ ' ' 

(B.7) 

and again using eq. ( |B.6|) for the asymptotic behavior, we obtain (ii) after straightforward 
algebra. 

Lemma B.3. Let A be the unique fugacity corresponding to given mean number 
density p G (0, l/2m), and Aqo = l/v^(so) = (1 — ■So)so ^™ with sq = (1 — 2mp)/[l — {2m — 
l)p]. Then for l/M = o(l), 

(i) p{l; A, M) = p^{l; Aoo) [1 + 0(//M)] 

(ii) Psih,... ,/. = /; A, M) = p,,oo(/i,... = /;Aoo)[l + 0(//M)] 
Proof. According to eq. ( p.lO|) 

p(/;A,M) = Afi(A,/-l)^^^^^^^ (B.8) 

where Q{X, M) = Yf^o Cj(A)sj(A)-^ (see Lemma B.l and the discussion after eq. ( p.ll|) ). 
In order to find the asymptotic behavior for M — > oo, one has to keep in mind that A 
depends on M. Let us denote by Am the fugacity corresponding to the mean occu- 
pation probability p in a system of (finite) size M, i.e. the unique solution oi p = 
\M~^d\ogVL{\, M) / dX (see Corollary A.l). In leading order Am approaches Aoo as 



\m — Xr 



(B.9) 



where 6 is a constant independent of M. The proof of this relation can be worked out 
by writing Am = Aoo + ca/ with limM^oo ca/ = 0,| and inserting this into the determining 
equation for Am- Careful expansion of the coefficients Cj(AM) and the powers Sj(AM) 



-M 



■^Note that, according to Corollary A.l we have for Aoo — l/¥'(so), liniAf^oo -^(Aoo, A'/)/Af = p, i.e. 
because of the one-to-one correspondence between p and A it must hold limAf^oo Am = Aoo- 
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with respect to ejv/ in the expression for Q{\m,M) then yields limM^oo^MM = bXoo and 
lirriM^ooi^M — bXoo/M]M'^ = const.. 
Using eq. ( p.9|) we find 



Q{Xm, M -l)= exp{bp) co(Aoo) So(Ac 



1 + 



M 



(B.IO) 



and thus obtain fi(A, M - l)/n{X, M) = So{X^y[l + 0(1 /M)]. Since Xm^{Xm: I ~ 1) — 
Aoo^^(Aoo,/ - 1)[1 + 0(//M)], it follows (i) from eq. (p^). Analogously, (ii) follows by 
using eq. ( ^.91 ) and the asymptotic expansion ( |B.1(J| ). 

Theorem B.l. If (M — 1) — (4m — 2)N > then the occupation probability 
p{l; N, M) is independent of / for all / G 7^l = {/ G N | /^^^ < Z < M + 1 - Z^^)}, with 
Z(^) = (2m - 1)A^ + 1, i.e. we can write p{l;N,M) = u{N, M) /Z{N, M). At the outer 
boundary points = l^^^ - 1 and = M + 2- p{l^^'> ; N , M) is different from 
p(/; iV, M) inside 7^l, in particular p{l^^^;N, M) = [u{N, M) + {-l)^]/Z{N, M). 

Proof. Given p(Z; A^, M) from eq. ( |2.6|) we first show that 

^ ^/ - (2m - l)(r + 1)' 



// 

m(/; at, M) = ^ Z{r -1,1- l)Z{N -r,M-l), Z{r, I) = i 

r=l ^ 



:b.ii] 



is constant inside TZi as long as (M — 1) — (4m — 2)A^ > 0. To this end we will proof that 
u{l]N, M) for / G TZi can be rewritten by use of the following combinatorial identity 



N 



u{l;N, M) 



E 



r=l 



l-l-{2m-l)r\ [ M -l-{2m-l){N +l~r\ 

N-r 



r — 1 

M+l-4m-(2m-l)(A^-l 
A^-1 

N 

-EM 



r=2 

= u{N,M) 



2m(r-l)-l\ /M+l-4m-(2m-l)(A^-r; 



r — 1 



N-r 



(B.12) 



which is independent of /. 

To verify the combinatorial formula we follow the methods described in the book of 
Riordan on combinatorial identities and use the following theorem of Lagrange for 
implicit functions [^: Let (f){z) be a power series in z with 0(0) 7^ 0, and let z{t) be the 
unique power series with z{0) = satisfying the implicit equation z{t) = t(j){z{t)). Then, 
for any power series F{z)., 



Fiz) 



'I - Wiz)) ^ n\ dX 

^ " n=0 



J-Tl jJTl 

E--LmA)$(A)niA=o 



Applying this theorem to (f){z) 
one obtains 



[l + z)-/^ and F( 



(B.13) 

1 + z)" with 1^1 < 1 and G N 



[1 + z' 



l + {(3 + l)z 



^^a;(a,n)t", u!{a,n) 

n=0 



< n < [f ] 

/3J 



(B.14) 
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Using (|B.14D we can write 



,0+7+2 



y + z) 



l + {(3 + l)z 



[l + z 



|a+7+l 



^ cj(0, n)r ^ + 7, n)e 



n=0 



n=0 



:b.i5) 



and 



(1 + ^) 



Q+7+2 



a+l 



:i + z 



|7+1 



X^a;(a,n)rX;^(7,^)i"- (B.16) 



n=0 



n=0 



_1 + (/?+ 

Comparing (p.l5| ) with (p.l6| ) and equating expansion coefficients we obtain 

n n 

^^0^(0, k)tu{a + 'j,n — k) = uj{a, k)uj{'j, n — k) . 



(B.17) 



fc=0 



fc=0 



In particular, for n < min{[a//3]; [7//^]} (see eq. (|B.14|) ),P| 



E 

k=0 



a — pk\ (1 — Pin — k) 
n — k 



k 

a + 'J — j3n 
n 



By setting n = N — 1, k = r — 1, j3 = 2m — 1, a = / — 2m, and 7 = M + 1 — 2171 — I for / G 7?.i 
(i.e. for (2m-l)A^+l < / < M-{2m-l)N), we have A^-l < min{[|^], [ ^Vm-^'' ]}' 
and can use eq. ( |B.18| ) to get eq. ( |B.12| ). 

In order to show that TZi is a "maximal" set, we again use eq. ( [B.17| ) with a = l^'^^—2m, 
/? = 2m-l, 7 = M + l-2m- Z^^), n = A^-landr = fc-l (note that n-l< [a/P], 
n > [a//3], n < [7//?], uj{M + 1 - 2m - 0) = 1 and uj{l^+^ - 2m, N - 1) = (-1)^"^) 
to obtain 



M(/W;iV,M) =u{N, M) + 



\N 



:b.i9) 



This completes the proof of theorem B.l. 

Corollary B.l. For M + 1 > 6m (that means there can be more than just one rod 
in the system) there does not exist a central region where p{l\ A, M) is constant (except 
for the trivial set Hi = {M/2, M/2 + 1} for even M). 

Proof. Assume that there exist a non-trivial central region TZi = {Iq, ... , M + 1 — Zq} 
with 2m < /o < [M — l)/2 in which p{l] A, M) is constant. Then there exists a function 



/(A, M) = Afi(A, / - l)n(A, M-l) 
independent of / for all I ^TZi. 



(B.20) 



^Equation ([B.18 ) may be identified as a combinatorial identity derived explicitly in the book of Riordan 
on page 148 ||27| by replacing f3 by — /3. However, in order to do this one has to define the binomial 
coefRcients for negative entries by analytical continuations of the Gamma function. 
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From 

Xn{X, I - l)fi(A, M-l) = Y^ Z{N, M)p{l, N, M)\^ , (B.21) 

N=l 

for A^o = max{iV > 1 : Z{N,M) > 0} = [(M+l)/2m]-l (note that Z{N,M) = (M-{2m-i)(7v+i)^^^ 
we get 

for all / G TZi and all 1 < < A^^q. But, according to theorem B.l, we have p{l, N, M) = 
const, in a central region if and only if / G TZi. Therefore, 

TZi C TZi = {{2m. - 1)N + 1, . . . , M - (2m - 1)A^} for all 1 < N < No . (B.23) 

Choosing N = Nq (and because we require TZi to have more than two elements for 
excluding trivial situations), we obtain 

(2m - l)No + 4 < M + 1 - (2m - l)No , (B.24) 

from which one readily concludes that M + 1 < 6m — 1 in contradiction to the restriction 
imposed on M. 

Comment. One can even show [|l^] that there does not exist any non-trivial region 
inside which p{l; A, M) is constant. 

Theorem B. 2. If M-l-(4m-2)A^+2(m-l)(s-l) > then the joint probabilities 
Ps{li, . . . ,ls',N,M) are translationally invariant for {h, . . . ,ls) G T^.^ = {(/i,... ,1^) G 
N^|/W < h; 2m < h-h-i hr k = 2,... ,s;ls< M+l-l^'^} with = (2m-l)(A/+l-s)+l, 
i.e. there exists a function f{yi, . . . , y^-i; N, M) exhibiting the property 

p,(/i, . . . , iV, M) = f{h - /i, - /.^i; AT, M) . (B.25) 

If (/i, . . . , /s) ^ Tig and (/i, . . . , Zt- ± 1, • • • , ^s) G 7?.s for some r G {1, . . . , s} then 

p,(/i, AT, M) = /(/2 - /i, - /.-i; AT, M) + {-lf+^-'/Z{N, M) . (B.26) 

Proof. The joint probabilities j9(/i, ... ,ls', N, M) are given by (see eq. ( |2.8| )) 

Ps{h,... ,ls;N,M) = ^ [[Z{rk-rk-i-l,lk-lk~i-l) 

l<ri<...<rs<N ^ ' ' k=2 

(B.27) 

By introducing new variables xi = ri, Xk = Vk — rk-i for A; = 2, . . . , s (X]fe=i = rg) we 
can rewrite this as 

n J .r ... Z{x,-l,h-l)Z{N-j:i,x„M-h) ^ ^, 

= Z(Ar,M)-i 5^ I f[Z(xfc-l,/fc-4_i-l) X 

X ^ Z(xi-l,/i-l)Z(A^-^a;i-xi,M-Z,) I , (B.28) 

2;i=l i=2 J 

20 



where An,s = {{xi, . . . , x^) G N*) | Xi + . . . + < A^}. 

Setting N' = N — J2i=2^i = ^ ~ {h — h), we have from theorem B.l (see 

eq. (TO ) 



AT' 



J2 Z{xi-l,h-l)Z{N'-xi,M'-h) = u{h;N',M') , (B.29) 

Xl = l 

where u{h;N',M') = u{N',M') = const, for {2m-l)N' + l <h< M' - {2m-l)N' . For 
the Xi this means (2m-l)(Ar-X]-=2 < < M-(/,-/i)-(2m-l)(A^-^'^2 ^i). The 

latter inequahty holds true for all (x2, • • . , x^) G ^Ar_i,s_i, if (2m — 1)(A^ + 1 — s) + 1 < li 
and/, < M-(2m-l)(iV+l-s), i.e. for(/i,... G 7^,. Since (Is-h) = Yfk=Sk-lk-i) 
we obtain 

= /(/2-/1,... ,/.-/s-i;iV,M), (B.30) 

which completes the proof of the first part of the theorem. 

To prove the second part, we insert eq. ( p.29|) in eq. (|B.28|) and consider an outer 
boundary point with (/i, . . . , ^ TZs and (/i, ...,/,- ± 1, G T?,, for some r G 
{1, . . . , s}. In fact, according to the constraints implied by the finite rod lengths, one can 
show that only r = 1 and t = s are possible. Due to symmetry we can restrict ourselves 
to the case r = 1. Then the outer boundary point is (Zi = l^^^ — l,/2; • • • ih) and we 
obtain 

p,(/i = /W-l,... ,/,;iV,M) = ^ \\{Z{xk-l^k-lk-i-l)u{luN\M' 

^ ' <Ajv-i,s-iU=2 ; h=lM-l 

(B.31) 

Except for the particular configuration (x2 = 1, . . . , x, = 1), all configurations (x2, . . . , x,) G 
•An-i,s-i yield arguments {li,N',M') for which u{li, N' , M') is constant (see the discus- 
sion above). For (x2 = 1,... ,Xs = 1), h is an outer boundary point of the set TZi 
corresponding to a system of size M' with N' rods. According to theorem B.l we thus 
have u{li, N', M') = u{N', M') + (—1)^ . By inserting these results in eq. (|B.31|) and by 
using the definition of /(yi, . . . , Us-i] N, M) in eq. ( |B.30|) we obtain 

p,(/i = /W -1, iV, M) = f{h - /(^) + l, Is-u N, M) + ^^^^ . (B.32) 

Corollary B.2. For (/i, ...,/,) ^ = {{h, ...,/,) G | 2m < /i < 4m, 2m < 
h—h-i for A; = 1, . . . , s; M+1— 4m < < M+1— 2m} (see footnote 1 in Sect. |3.2| ) there 
does not exist a region, where Ps{h, • • • ,h] ^, M) is translationally invariant. 

Proof. Ifpsih,... , Is] A, M) = f{l2 -h,... Js- Is-i] A, M) for some {h, . . . ,ls) ^ Cs 
then, from eq. ( p. 9]) , 

n(A, /i - 1)(](A, M' - /i) = ^^i^2-^/(/2 - /i, - A, M) 



(B.33) 



21 



where M' = M — (Ig — h). Accordingly, there exists a range of consecutive li values, 
where fi(A, li — l)i7(A, M' — li) is independent of li. This, however, is impossible due to 
the comment after Corollary B.l. 

Theorem B.3. If M - 2(2m - 1)(A^ + 1 - s) > then the joint probabilities 
. . . , Is] N, M) are constant functions for (/i, . . . , Z^) G i3s = {(/i, . . . , G N'^ | (2m— 
l)(iV + 1 - s) + 1 < /i ; (2m - 1)(A^ + 1 - s) + 1 < /fc - h-i for = 2, . . . , s ; < M - 
(2m-l)(A^ + l-s)}, i.e. we can write p,(/i, .. . ,ls;N,M) =v{N,M). li {h, . . . ,/,) ^ 
and (Zi, ... , ± 1, . . . Js) G i3s for some r G {1, . . . , s} then 

p,(/i, iV, M) = v{N, M) + ^^^^ . (B.34) 

Proof. The explicit formula (see eqs. 



Z{N,M)ps{h,...,h;N,M) = C' ^^''') >< (B-35) 

(xi,... ,Xs)£Am,3 

M+l-2m-h-{2m-l){N yr //fc-4-i-l-(2m-l)(a;fc-a;,.__i; 

can be rewritten by using eq. ( p.l4|) derived in theorem B.l (remember that we defined 
An^s = {{xi, ... ,Xs) eW) I Xi + . . . + < A^} after eq. ( PT^ )). We have for > 0, 
i = 1, . . . , s + 1, and £n,s = {("-i, • ■ ■ , Ug) e Nq \ rii + . . . + Us = n} , 

^ fc=lnfc=0 n=0 (ni,...,ns+i)G£-„,s+i fc=l 

(B.36) 

and 

nrrkTHMrVk^ ^5-<- + ■ ■ ■ + — n E 

oo s+1 

= J] Ecc;(«i + . . . + ni) Y[ uj{0, Uk) . (B.37) 

n=0 fn.s+i fe=2 

Comparing ( |B.36| ) and ( |B.37| ) and equating coefficients we obtain 

s+1 s+1 

^ JJu;(afe,nfc)= ^ t^(ai + . . . + a^+i, ni) JJ u;(0, n^) , (B.38) 
which, for n < mini<fc<s+i{[a;fc//3]} (see eq. ( p.l4| )), yields 

We have defined (^^) = 1 on the right side here (but not on the left hand side). 
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By choosing ai = h — 2m, ak = Ik — h-i — 2m for A; = 2, . . . , s, Ug+i = M + 1 — — 2m, 
j3 = 2m — 1 and n = N — s < mini<k<s+i{(^k/ P} , we can identify the combinatorial 
expression (p.36| ) by the left hand side of (p.39| ). Substituting the right hand side we 
then obtain 

p.(/„...,/.;iv,M) = ziN,MrYl (-ir^-^^'P'^^lS^^'^^^^V 

'M-2ms-(2r 



2m-l)xi\ -Q(-_^^|a.,-i/^2m(xfc-l)-l\ 



^ ^ fc=2 

= v{N,M), (B.40) 

which completes the proof of the first part of the theorem. 

Now, let (/i, ... , /s) ^ -Bs and (/i, . . . , ± 1, . . . , Z^) G i3s for a given r G {1, . . . , s}, 
and /o = and Z^+i = M + 1. Then - s < [afc/(2m - 1)] for A; 7^ r and iV - s - 1 = 
[a^/ (2m — 1)], which can be shown by simple topological considerations. Accordingly, for 
(ni, . . . , n^+i) G = ^Af-s,s+i\('^i = 0, . . . , = iV-s, . . . , n^+i = 0), ci;(afc, n^) = 

(°''~*-^™~^-'"''"') , while for the particular element {n\ = 0, . . . ,nr = N — s, . . . , Ug+i = 0), 

iu{ak,nk) = (1 — Skr) + {—^Y^~^^SkT- Evaluating the left hand side of eq. p.38| we arrive 
at 

E fi^M = E n("-p';^-i'M + (-i)- 

(ni,... ,ns+i)G£']v-s,s+i !^=^ (ni,... ,ns+i)e£'j^_^ ,.^-^ k=l 

= Z{N,M)p,{h,... ,h;N,M) + {-lf-\ (B.41) 

On the other hand, the right hand side of eq. ( |B.3^ ) is equal to Z{N, M)v{N, M) and 
hence we obtain 

p,{h, ...,h-N,M)= v{N, M) + ^-^^j^ ■ (B.42) 



APPENDIX C 

For < M' + a < M and a > 0, fi(A; M', a) is defined as follows (for v{n) = VuM, 
see Sect. |) 

fi(A; M', a) = 0, for M' < 2m - 1 

n{X;M',a) = 1, for 2m-l < M' <4m-l (C.l) 

A/'+l-2m 

fi(A;M',a) = 1 + A ^ e-"("+")n(A; r - 1, a), forM' > 4m - 1 

r=2m 

Theorem C.l. The occupation probability 



23 



is the unique solution of the following set of nonlinear coupled difference equations 

l+2m-l l+2m~l 

= -logX + u{l)+logPi+ log[l-t™(s)+p,]- log[l-tm(s)], (C.3) 



s=l+l 



s=l 



for / G {2m, . . . , M+1 — 2m} and pi = else, where im{s) = Yl'j=o ^ Ps-j- 

Lemma C.l. Let M' > 2m. If < M' + a < M and a > 0, then fi(A; M', a) obeys 
the recursion relations 

(i) Q{X; M', a) = Q{X; M' -l,a) + A-"(°+^^'-2™+i)r](A; M' - 2m, a), 

(ii) fi(A; M', a) = n{X; M' - 1, a + 1) + Ae-"("+2m)^(^. ^> _2m,a + 2m). 

Lemma C.2. For s > 2m let 

fi(A;s- 1,0) 



1 trn. i S] 



-<P{X-M,s) 



n{X;M, 0) 

Then 0(A; M, s) obeys for ps = p{s; X, M) (from eq. ( p.2|) ) the relations 

(i) 0(A; M, s) = n{X; M - s + 2m - 1, s - 2m + 1), 

(ii) 0(A; M,s-1) = 0(A; M, s) + Ae-"(")^](A; M -s,s). 
Proof of Lemma C. 1 . For M' > Am — 1 

M'-2m 

Vt{X-M',a) = 1 + A ^ e-"("+")^](A;r-l,a) + Ae-''("+^'-^"^+^)^](A;M'-2m, 



(C.4) 



a 



r=2m 



n{X;M' - l,a) + e 



~u{a+M'-2m+ 



i)fi(A;M' - 2m, 



a 



(C.5) 



If 2m < M' < Am - 1, then n{X; M' -l,a) = 1 and n{X; M' - 2m, a) = 0. Hence, 
eq. (|C5|) is also valid for 2m < M' < Am — 1, which implies (i). 
According to (^), 



Z{N, M',a) 



N 



N 



X^exp- 

l<il,...<i]v<M' 
M' 



y^^u{a+ik) + v{ii) + y^^v{ik-ik-i)+v{M' + l-iM) 



k=l 



k=2 



1=1 



l+l<i2,...<iN<M' 



N 



N 



M 



1=1 



k=2 



-V U X 



fc=3 



7V-1 



Yu{a+ik) + v{i2-l) + y^^v{ik-ik-i)+v{M' + l-iN) 

k: 

exp — [u{a + l) + 

N-: 

^u(a + /+jfc)+t;(ji) + ^t;(jfc-jfc_i) + i;(M'-/+l-zjv) 

(C.6) 



k=l 



l<ii<-<j]V-i<A/'-i 
M' 



k=2 



1=1 
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Remembering that Z(0, M' , a) = e (^see the remark right after eq. ( |2.3| )), it holds 

oo 

Q{X,M',a) = 

N=0 

M' 

= e-"(^'''+^) + A ^e-"("+^)e-"(')fi(A,M' -/,« + /), (C.7) 



1=1 



where v{n) = oo for n < 2m, and v{n) = for n > 2m. Since fl{X, M' — l,a + l) =0 for 
M' — I < 2m — 1 we have 



M' 



n{X, M', a) = 1 + A ^ e-"("+')n(A, M'-l,a + l), (C. 



l=2m 

which for r = M' — I yields 

M'-2m 



1 + A J2 e-"("+^'^'-'')l](A; r, a + M' - r) . (C.9) 



r=2m—l 

Hence we finally obtain (ii), 

M'-2m 



Q{X,M',a) = 1 + A ^ e-"("+^^'-'')(](A;r,a + M'-r) 



r=2m—l 
M'-l-2m 

= 1 + X e-"«"+iW^'^'-i>-)fi(A;r,(a + l) + (M'-l)-r) 

r=2m—l 

+Ae-"(°+2'")fi(A; M'-2m, a + 2m) 
= fi(A,M'-l,a + l) + Ae-"("+^")l](A;M'-2m,a + 2m). (C.IO) 

Proof of Lemma C.2. We prove the proposition (ijhj complete induction with respect 
to s. For s = 2m we have 

^(A;M.2,„) = ^^'^|^;^^°)^^ ll-U2.n)]. (C.ll) 

With fi(A; 2m — 1, 0) = 1, im{2m) = Ylj=o ^ P2m-j = p2m (note that pi = for i < 2m, 
and = Ae-"(2™)J1(A; M - 2m, 2m)/r](A; M, 0)) it follows 

0(A; M, 2m) = ^](A; M, 0) - Ae-"(2™)^](A; M - 2m, 2m) . (C.12) 

For M' = M and a = we then obtain by using Lemma C.l. (ii) 

Q{X, M, 0) = Q{X, M - 1, 1) + Ae~"(2'")fi(A; M - 2m, 2m) (C.13) 

and hence 

0(A; M, 2m) = r](A, M, 0) - Ae-"(2"')(](A; M - 2m, 2m) = fi(A, M - 1, 1) . (C.14) 



25 



Accordingly, proposition (i) is valid for s = 2m. Let us now assume that (i) holds true 
for s — 1 > 2m — 1. 

Since im{s) = E5=o"^P«-j and pi = Ae-"WQ(A;Z - 1,0)Q(A;M - Z)/Q(A; M, 0) we 
can write 



<P{X;M,s) 



Q{\;M,0) 





1,0) 


1 




0(A;s- 


1,0) 


1 





-u(^) MX; s-J-l, OMX; M-s+j,s- 
n(A;M,0) 

2rTV-l 

(}(A;M, 0)-A J2 e~''^^^n{X;s-j-l,0)n{X;M-s+j,s-j) 

j=o 

2m-l 

Q{X] M,0)-XJ2 e-'''^'^^~^^n{X] s-j-2, 0)Q{X] M-s+l+j, s-1-j) 



Q(A;s-l,0) 
A 



j=0 



+ 



[e""^"~'^""^^{X; s-2m-l,Q)n{X; M -s + 2m, s-2m) 
ilyX] s — 1, Oj 

-e-"(")Q(A; s-1, 0)Q(A; M-s, s)] 
n{X;s-2,0)<j){X;M,s-l) 



n{X;s-l,0) 

+ Ae-"(^2"*)n(A; s-2m-l, 0)Q(A; M-s+2m, s-2m)l -Ae-"(^)Q(A; M-s, s) 

where the line 3 follows from line 2 after some straightforward manipulation of the sum 
over j, and line 4 from line 3 by definition of 0(A; M, s — 1). 

Due to the induction hypothesis, 0(A; M, s—1) — Q{X; M—s+2m, s— 2m), from which 
follows 



(C.15) 



0(A;M,5) 



fi(A:.l/-,s + 2///.,s-2m) 



O(A;s-l,0) 
-Ae-"(^)Q(A;M-s,s). 



[n{X; s-2, 0) + Ae-"(^^'")n(A; s-2m-l, 0)] 

(C.16) 



By using (i) of Lemma C.l. for M' = s — 1 and a = we obtain 

Q(A; s - 1, 0) = n(A; s - 2, 0) + Ae-"(^-2"*)n(A; s - 2m - 1, 0) (C.17) 

and hence 

(/)(A; M, s) = 0(A; M - s + 2m, s - 2m) - Ae-"(^)0(A; M - s, s) (C.18) 

By using (ii) of Lemma C.l. for M' = M — s + 2m and a — s — 2m we find 

n{X; M -s + 2m,s-2m) = n{X; M - s + 2m - 1, s - 2m + 1) + Ae-"(")f)(A; M-s. s) 

(C.19) 

and hence 

Q(A; M - s + 2m, s - 2m) - Ae~"(^)Q(A; M - s, s) = Q(A; M - s + 2m - 1, s - 2m + 1) 

= (f){X;M,s). (C.20) 

This completes the proof of part (i) of Lemma C.2. 
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To proof the second part (ii), we apply (ii) of Lemma C.l. For M' = M — s + 2m and 
a = s — 2m we obtain 

n{\; M - s + 2m,s -2m) = Q{X; M - s + 2m + 1, s - 2m + 1) + Ae-"(")^](A; M -s,s). 

(C.21) 

Due to part (i) of the Lemma just proven, 0(A; M, s) = f2(A; M—s+2m — l, s— 2m+l) and 
0(A; M, s — 1) = Q{X; M — s + 2m, s — 2m). Inserting this in ( |U.21| ) it follows proposition 
(u). 

Proof of Theorem C. 1 . 

According to Lemma C.2. and the definitions of Ps and tm{s), 



1 - tm{s) +Ps 



Q{X;s- 1,0) 
n{X;M, 0) 

n{X;s- 1,0) 
n{X;M, 0) 



[0(A; M, s) + Ae-"(^)fi(A; M - s, s)] 
0(A;M,s-l). (C.22) 



Since fi(A; s - 1, 0) = for s < 2m and fi(A; M - s, s) = for s > M + 1 - 2m, it follows 
ps = for s ^ {2m, ... , M + 1 — 2m}. For s G {2m, ... , M + 1 — 2m} on the other hand, 
we obtain from Lemma C.2. and eqs. ( p.2| , |C.3| , |C.22| ) 



= ^ogA+M(/)+log 



s=l+l 



^^_,(,)n(A;s-l,O)fi(A;M-/,0 
n{X;s-l,0) 



fi(A;M,0) 



-0(A;M,s-l) 



l+2m-l 



fi(A;M,0) 

log fi(A; M -1,1)- log 0(A; M, / + 2m - 1) 
0. 



s=l 



fi(A;s-l,0) 
_ f](A;M,0) 



0(A;M,s) 



(C.23) 
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Figure Captions 



Fig. 1 Occupation probability Poo{l) of hard rod centers as a function of the distance / from 
a hard wall for (a) three small rod lengths and a large mean occupation number 
p = 0.1, and (b) three large rod lengths and a small mean occupation number 
p = 0.02 corresponding to a continuum- like situation. The solid lines in (a) were 
drawn as a guide for the eye. 

Fig. 2 Correlation function C{1) between the first possible position 2m of a rod center and 
another rod center that is at distance / from the wall (see eq. |4.2| ) )for the same 
parameters as in Figs. la,b. The solid lines in (a) were drawn as a guide for the eye. 

Fig. 3 Occupation probability Poo{l) of rod centers as a function of the distance / from a 
wall for (i) hard rods in the presence of a hard wall {v = Vhr, u = 0), (ii) hard rods 
in the presence of a soft wall {v = Whr, u = Uq), and (in) rods with a Lennard- Jones 
type Takahashi interaction in the presence of a hard wall {v = Wlj, u = 0). In (a) 
the discrete nature of the lattice is important {m = 4, p = 0.1), while in (b) the 
data correspond to a continuum-like situation (m = 18, p = 0.02). The solid lines 
in (a) were drawn as a guide for the eye. 

Fig. 4 Correlation function C{1) between the first possible position 2m of a rod center and 



another rod center that is at distance / from the wall (see eq. 4.2) for the same 



parameters as in Figs. 3a,b. The solid lines in (a) were drawn as a guide for the eye. 
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FIG. 1. Occupation probability Poo(0 hard rod centers as a function of the distance 
I from a hard wall for (a) three small rod lengths and a large mean occupation number 
p = 0.1, and (b) three large rod lengths and a small mean occupation number p = 0.02 
corresponding to a continuum-like situation. The solid lines in (a) were drawn as a guide 
for the eye. 



(Fig.l, Buschle et al.) 



31 



0.08 

Cf/j 0.04 
0.00 
-0.04 
-0.08 
-0.12 

0.02 

C(l) (b) 

0.01 ^ ' 

0.00 
-0.01 
-0.02 




10 20 30 40 50 60 



p=0.02 





100 



200 300 

1/ 



m=14 
m=18 
m=22 



400 500 



FIG. 2. Correlation function C{1) between the first possible position 2m of a rod center 
and another rod center that is at distance / from the wall (see eq. |4.2| ) )for the same 
parameters as in Figs. la,b. The solid lines in (a) were drawn as a guide for the eye. 



(Fig. 2, Buschle et al.) 
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FIG. 3. Occupation probability Poo(0 of rod centers as a function of the distance / from 
a wall for (i) hard rods in the presence of a hard wall {v = Whr, m = 0), (ii) hard rods in 
the presence of a soft wall {v = Whr, m = Uq), and ('ui^ rods with a Lennard- Jones type 
Takahashi interaction in the presence of a hard wall {v — Vi^j, u — 0). In (a) the discrete 
nature of the lattice is important {m — A, p = 0.1), while in (b) the data correspond to a 
continuum-like situation (m = 18, p = 0.02). The solid lines in (a) were drawn as a guide 
for the eye. 



(Fig. 3, Buschle et al.) 
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FIG. 4. Correlation function C{1) between the first possible position 2m of a rod center 



and another rod center that is at distance / from the wall (see eq. [4 .21) for the same 
parameters as in Figs. 3a,b. The solid lines in (a) were drawn as a guide for the eye. 



(Fig. 4, Buschle et al.) 
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